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Abstract. We calculate multicomponent stellar wind models with 
inclusion of a helium component applicable to He-rich and He-poor 
stars. We show that helium does not decouple from the stellar wind 
of He-rich stars due to its coupling to hydrogen and its ionized state. 
For He-poor stars helium may decouple from the stellar wind, how- 
ever this effect is not able to explain the chemical peculiarity of these 
stars. We conclude that the explanation of chemical peculiarity of 
these stars based purely on helium or hydrogen decoupling from the 
^ ; stellar wind is unlikely. 
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1 Introduction 

Radiative force may have important influence on the structure of stellar atmospheres. Al- 
though there is not known any significant direct influence of the radiative force on the atmo- 
spheres of most cool stars, the radiative force may be dominant in the atmospheres of hot 
stars. The consequences of this influence are different. Mainly for A stars and white dwarfs 
the radiative force may cause the radiative diffusion and subsequent elemental abundance 
anomalies (see Alecian 1995 or Vauclair 2003 for a review). On the other hand, the radiative 
force may accelerate a stellar wind mainly for O stars and WR stars (Kudritzki & Puis 2000). 
However, there exists a group of B stars, for which both effects of stellar wind and elemental 
diffusion in the stellar atmosphere are important. It is generally believed that in order to 
explain the elemental abundances of helium and some other elements in the atmospheres of 
He-rich stars at least three different ingredients are necessary. These ingredients are stellar 
wind, elemental diffusion and magnetic field (e.g. Michaud et al. 1987). 

However, alternative explanation of chemical peculiarity of He-rich and He-week stars 
exists. It was proposed by Hunger & Groote (1999, hereafter HG). They showed that for 
stars with effective temperatures T e s < 25 000 K helium may decouple from the stellar wind, 
fall back onto the stellar surface and create regions with enhanced abundance of helium if 



2 



magnetic fields are present. Similarly, for stars with effective temperatures T e g < 17 000 K 
even hydrogen may decouple from the stellar wind, fall back onto the stellar surface and cause 
helium underabundance if magnetic fields are present. We decided to test this explanation of 
chemical peculiarity of He-rich and He-week stars using our multicomponent wind models. 
First preliminary calculations have already been performed by Krticka & Kubat (2001a), 
however with artificially low helium charge and without ionization balance calculation. 



2 Multicomponent stellar wind models 

Stellar winds of hot stars are accelerated mainly by the absorption of radiation in the res- 
onance lines of such elements like carbon, nitrogen, oxygen, or iron. However, these wind 
components have much lower density than the rest of the stellar wind, which is composed 
mainly of hydrogen and helium. Thus, the process of acceleration of a stellar wind of hot 
stars has two different physical steps. In the first step, radiation is absorbed by low-density 
metals like carbon, nitrogen, oxygen, or it is scattered by free electrons, and these wind 
components obtain momentum from stellar radiation field. In the second step, momentum 
obtained by these low-density components is transferred to high-density wind components 
via the Coulomb collisions between charged particles. Thus, stellar winds of hot stars have a 
multicomponent nature. For stars with relatively high wind density (e.g. galactic O stars or 
WR stars) this multicomponent wind nature does not influence wind structure, as has been 
discussed by Castor et al. (1976). However, for stars with a relatively low wind density (i.e. 
main-sequence B stars or stars with extremely low metallicity, see Krticka et al. 2003) the 
multicomponent wind nature becomes important for the wind structure itself because new 
effects occur, for example frictional heating or decoupling of wind components (Springmann 
& Pauldrach 1992, Babel 1995, Krticka & Kubat 2001b, hereafter KKII). 



2.1 Model equations 

Equations used here for the calculation of four-component wind models are nearly the same 
as those of KKII. We assume stationary and spherically symmetric stellar wind, which is 
composed of four components, namely absorbing metallic ions, hydrogen, helium, and free 
electrons. For the calculation of wind models we solve the continuity equation, momentum 
equation, and energy equation for each component of the flow. The continuity equation has 
the form 

{r 2 paV ra ) = 0, (1) 

where p a is the density, r is the radius, and v ra is the velocity of a component a. The 
momentum equation is 

d^rq _ rad 1 d ( 2 n \ , la „ f ric /r>\ 

V ra^T-9a ~ 9 ~ J ^ M + — E + ^ 9ab , (2) 

where g™ A is the radiative force, E is electric polarisation field, and g^ b lc is the frictional force 
(Burgers 1969) 

fric 1 ts sit \ V ' rb ~~ Vra /o\ 

9ab = —K ab G{x ab )- (3) 

Pa \Vrb-V ra \ 
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where G(x a b) is the Chandrasekhar function, and the frictional coefficient is 

K ab = n a n b ^0l\nA, (4) 
kT ab 

where mean temperature T ab = (m b T a + m a T b ) / {m b + m a ) is calculated using temperatures 
T a and T b of individual wind components with atomic masses m a and m b . The radiative 
force in the CAK approximation (Castor et al. 1975, Friend & Abbott 1986, Pauldrach et 
al. 1986, see also KKII for generalization of CAK force for a multicomponent flow) due to 
line-absorption acts on metals and radiative force due to the Thomson scattering acts on free 
electrons. 

Energy equation for each component of the flow is 

\v ra ^ + aJ ^i (*0 = <T +£(«+Q£<) ■ (5) 

where the heat exchange is given by 

1 2k (T b - T a ) exp (-x 2 ab ) 

Qab = —^&ab : , (6) 

Vtt m a + m b a ab 

Qat = 7 K ab G(x ab )\v rb - v ra \, (7) 
m a + m b 

and Q r ^ d is the radiative heating calculated using the the method of thermal balance of 
electrons (Kubat et al. 1999). We also take into account the Gayley-Owocki (1994) heating. 
The base flux, which is necessary for the calculation of the radiative heating term, is taken 
from H-He spherically symmetric NLTE model atmospheres (Kubat 2003). 

System of hydrodynamic equations is closed using equation for electric polarisation field 
and by equations of ionization equilibrium (we assume nebular approximation after Mihalas, 
1978, Eq. 5.46). For a more detailed description of three-component variant of these models 
see KKII. 



frictional heating is 



2.2 The possibility of wind decoupling 



The frictional force Eq. (JH} between two components depends on the velocity difference be- 
tween these two components via the so-called Chandrasekhar function, 
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where non-dimensional velocity difference x ab is given by 
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where a ab is the mean thermal speed, a 2 ah = 2k {m a T b + m b T a )/ {m a m b ). The plot of the 
Chandrasekhar function is given in FigO For relatively low velocity differences, Av ab < a ab , 
the Chandrasekhar function is an increasing function of the velocity difference Av ab and 
components a and b are well coupled. However, for higher velocity differences, Av ab > a ab , the 
Chandrasekhar function is a decreasing function of the velocity difference and the decoupling 
of components may occur. 
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Figure 1: The run of Chandrasekhar function. Note that Eq. © yields x ab ~ Av ab , 
where Av ab is the velocity difference between wind components. If the flow is well coupled, 
Xab < 1-0, G(Av ab ) ~ Av ab . If the drift velocity is large, x ab > 1.0, G{Av ab ) ~ Av~ b 2 , and 
wind may decouple. Note that the point x ab ~ 1.0 corresponds to the maximum of G. 



2.3 Solution using Newton— Raphson method 

Multicomponent wind equations (I1I2I5D can be formally written as 

Pip = 0, 

where the vector describing the solution has the form of 



$ = (y>i,i/> 2 ,. . . ,v 
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where 
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The velocity difference Av r ^ was taken as an additional independent variable to improve 
convergence properties. The solution can be obtained using the iteration scheme 



where the Jacobi-matrix is 



rsif) n+1 = -p n v> n , 

rn _ dPk 



(13) 
(14) 



More details about the solution procedure can be found in Krticka (2003). 



3 Calculated models 



We have calculated four-component stellar wind models (i.e. models consisting of absorbing 
ions, hydrogen, helium, and free electrons) to test the explanation of helium chemically 
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Figure 2: Four- component wind model of a main-sequence B3 star. Upper left 
panel: Wind velocities of wind components. Note that velocities of all wind components are 
nearly the same. Upper right panel: Temperatures of wind components. Temperatures of all 
wind components are close to each other. Note slight heating due to the Gayley & Owocki 
(1994) and frictional heating. Lower panel: Charges of wind components. 



peculiar stars proposed by HG. Here we present two stellar wind models of B3 and B6 main- 
sequence stars. Spectral type of helium-rich stars is close to B3, for which HG predicted 
helium decoupling from the stellar wind. On the other hand, spectral type of helium-poor 
stars is close to B6. For these stars HG predicted hydrogen decoupling from the stellar wind. 

3.1 Models suitable for He-rich stars 

We have selected a main-sequence B3 star as a representative example of He-rich stars. Stellar 
parameters are taken according to Harmanec (1988) and force multipliers are after Abbott 
(1982). The plot of velocities, temperatures and charges of individual wind components is 
given in Fig.|21 

Contrary to the predictions of HG, helium does not decouple from the main wind. The 
reason for this behaviour is a neglect of helium coupling to hydrogen by HG. They assumed 
that helium is accelerated mainly due to the collisions with metallic (absorbing) ions. The 
frictional force fu e ,i = Pue9ne,i is proportional to (see Eqs. (|3l4j) ) 

,/He,i ~ ^He^i, 
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Figure 3: Four- component wind model of a main-sequence B6 star. Upper left 
panel: Wind velocities of wind components. Note that helium decouples from the stellar wind 
and subsequently decelerates. Helium velocity is lower that the escape velocity. Thus, helium 
may fall back onto the stellar surface. Upper right panel: Temperatures of wind components. 
Stellar wind is heated mainly due to the frictional heating in the outer parts of the wind. 
Lower panel: Charges of wind components. Note that helium recombines in the outer model 
region. 



where nu e and ri\ are number densities of helium and metallic ions. However, the frictional 
force due to the collisions between helium and hydrogen atoms is 

where % is the number density of hydrogen atoms. Since nu ^> n;, the frictional force due to 
hydrogen atoms may be important for the acceleration of helium atoms, and, consequently, 
can not be neglected. 

We conclude that helium does not decouple from the stellar wind of He-rich stars. Thus, 
the explanation of a helium overabundance in chemically peculiar stars by the helium decou- 
pling in their stellar wind is questionable. 



3.2 Models suitable for He-poor stars 

We have selected main-sequence B6 star as a representative for He-poor stars. Again, stellar 
parameters are taken from Harmanec (1988) and force multipliers are after Abbott (1982). 
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The plot of velocities, temperatures and charges of individual wind components is given in 
Fig.EI 

Apparently, in the outer model region helium atoms recombine. Because the frictional 
acceleration depends on the square of helium atomic charge (see Eqs. (jSEJ), the frictional 
acceleration is not able to support helium atoms any more and helium decouples from the 
flow. The decoupling is accompanied by large frictional heating. The helium velocity is lower 
than the escape velocity, thus helium may fall back onto the stellar surface and create helium 
overabundance in the stellar atmosphere. However, this spectral type corresponds to He-poor 
stars, thus even in this case the four-component models are not able to explain their chemical 
peculiarity. 



4 Conclusions and discussion 

We have calculated four-component wind models (i.e. models with absorbing ions, hydro- 
gen, helium, and free electrons) applicable to helium chemically peculiar stars. We used our 
models to test the explanation of helium chemically peculiar stars by helium and hydrogen 
decoupling in the stellar winds of these stars. For hot B stars (in the parameter range of 
He-rich stars) helium is coupled not only to metals but also to hydrogen atoms, and it is 
ionized. Consequently, helium does not decouple from the stellar wind of He-rich stars. We 
conclude that the explanation of chemical peculiarity of He-rich stars (HG) based on helium 
decoupling is questionable. 

On the other hand, for cool B stars (in the parameter range of He-poor stars) helium 
may recombine and consequently fall back onto the stellar surface. However, this decou- 
pling can not explain chemical peculiarity of He-poor stars, since it may produce a surface 
overabundance of helium, not an underabundance, as anticipated by HG. 

However, there are two effects which are still unclear. First of all, the mass-loss rates 
of B stars are highly uncertain. Our test calculations showed that due to lower mass-loss 
rate both hydrogen and helium may fall back onto the stellar surface in the case of He-rich 
stars. However, this effect can not help to explain the chemical peculiarity of these stars since 
hydrogen and helium are well coupled even in this case. Another problem is the calculation of 
wind ionization structure. From the theoretical point of view, as shown by Krticka & Kubat 
(2001), the helium decoupling is possible for artificially lowered helium charge. However, our 
simplified ionization calculations (based on a nebular approximation) do not predict such low 
helium charge. We think that our calculated ionization structure is roughly correct, but we 
shall perform tests using more detailed calculations. Both these effects, i.e. correct mass-loss 
rates and ionization structure, will be incorporated in our NLTE wind code. First results 
obtained using our code are currently under way (Krticka & Kubat 2003). 

We point out that helium diffusion superimposed with the multicomponent stellar wind 
may explain chemical peculiarity of He-rich stars (Michaud et al. 1987). On the other hand, 
mechanism of launching of He-free wind proposed by HG may possibly work for some cooler 
stars, however probably not for He-rich stars (see also Kubat & Krticka 2004). 
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